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Abstract. We have measured the angular correlation function, w{9), of radio sources in the 1.4 GIfz NVSS and 
FIRST radio surveys. Below ~ 6' the signal is dominated by the size distribution of classical double radio galaxies, 
an effect underestimated in previous studies. We model the physical size distribution of FRII radio galaxies to 
^\ ' account for this excess signal in w{6). The amplitude of the true cosmological clustering of radio sources is roughly 

constant at yl ~ 1 x 10""^ for flux limits of 3-40 mjy, but has increased to A ~ 7 x 10"'^ at 200 mjy. This can 
be explained if powerful (FRII) radio galaxies probe significantly more massive structures compared to radio 
^ I galaxies of average power at z ~ 1. This is consistent with powerful high-redshift radio galaxies generally having 

massive (forming) elliptical hosts in rich (proto-)cluster environments. For FRIIs we derive a spatial (comoving) 
correlation length of ro = 14 ± 3 Mpc. This is remarkably close to that measured fo r extre mely red objects 
I , (EROs) associated with a population of old elliptical galaxies at z ~ 1 bv iDaddi et al.l i200 j) . Based on their 

■ similar clustering properties, we propose that EROs and powerful radio galaxies may be the same systems seen at 

' different evolutionary stages. Their ro is ~ 2x higher than that of QSOs at a similar redshift, and comparable to 

, that of bright ellipticals locally. This suggests that ro (comoving) of these galaxies has changed little from z ~ 1 

to z = 0, in agreement with current ACDM hierarchical merging models for the clustering evolution of massive 
early- type galaxies. Alternatively, the clustering of radio galaxies can be explained by the galaxy conservation 
model. This then implies that radio galaxies of average power are the progenitors of the local field population of 
Q ' early-types, while the most powerful radio galaxies will evolve into a present-day population with ro comparable 

^ ! to that of local rich clusters. 

^ ' Key words, cosmology: large-scale structure of universe - galaxies: active - galaxies: statistics - radio continuum: 

galaxies - Surveys 
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1. Introduction dard for studying large-scale structure. The clustering of 

1 1 • 1 , 1 r cosmological objects can be characterized by their spatial 
In striking contrast with the extremely high level or , ,. r ,. i • i i ii r ^/ \ c / ^--v 

° ° correlation lunction, which has the lorm ^(r) = [r/ro) ' 

isotropy observed in the temp erature of the cosm icmi- ^^^^^ present-day correlation length and 7 1.8 

crowave background (see e.g. | de Bernardis et al. ||2000D, fo^ objects ranging from clusters to normal galaxies (see 

galaxies are not distributed throughout the Universe in a |Bahca'll fc Sone^rlfllQSa for a review). The local popu- 

random manner. According to the gravitational theory of j^^.^^ ^^^^^.^^ ^ relatively unbiased tracer of the 

instability the present structures originated from tiny flue- j 1 • 4.4. j- 4. -i, j.- -4.1, r a 1.-1 

, 7^11,1 underlying matter distribution, with tq = 5.4 n Mpc 

tuations m the initial mass density field. This has shaped derived from gal axies in the early CfA redshift survey by 

the large-scale structure of the Universe, which consists of Ip^yis fc Peeblej £983), however more recent low-redshift 

vast empty regions (Ws), and strings of dark and lumi- ^^^^^^^ ^^^^ ^^^^ clustering of galaxies depends 

nous matter (wa/Zs) where billions of galaxies are found. , , ji 1,1-14- i? 

^ ' strongly on luminosity and/ or morphological type. 1* or ex- 

The clustering properties of galaxies can be quantified ^^^^^^ ^^^^^ > elhpticals represent spatial structures 
using statistical techniques, such as methods of nearest ^^^^ ^^^^ more strongly cluste red with m - 7 - 12 
neighbours, counts m cells power spectra, and correlation j^^i |(^^^^^^^ IWillmer et alJ Hooi 
functions (see| PeebM | 1980t for an in-depth mathematical iNorberg et alJi2n02^ . From deep, magnitude-limited red- 
review). In particular the two-pomt correlation function ^^^^^ ^^^^p^gg ^^^^ ^^^^^ ^^^^ ^^^^ comoving corre- 
is a simple, but powerful tool that has become a stan- ^^^-^^ ^^^^^ galaxies declines with redshift, roughly as 
Send offprint requests to: R. A. Overzier expected from simple gravitational theory (e.g. CFRS, Le 
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Fevre et al. 1996; Hawaii K, Carlberg et al. 1997; CN0C2, 
Carlberg et al. 2000; CFDF, McCracken et al. 2001). In 
contrast to this, the clustering stre ngth of quasars ap- 
pears to vary little over < z < 2.5. ICroom et al.l l)200l[l 
found an approximately constant amplitude of ~ 5 
Mpc from ~ 10 OOP quasars in the 2dF QSO Redshift 
Survey. Likewise . IPaddi et all l)200lll2002(l found that the 
(comoving) correlation length of massive elliptical galax- 
ies also shows little evolution with redshift. They find 
To = 12 ± 3 Mpc for a popul ation of extremely red 
objects (EROs) at z ^ 1 (see also iMcCarthv et al ] l200lt 
iRoche et all l2002t iFirth et alJ l2002|) . which are consis- 
tent with being the pas sively evolving progenitors of local 
mass i ve ellipticals Ce.g. IPunlop et al.lll996t Cimatti et al.l 
ll998il2n02tlDev et al.lll999l:lLiu et al."2000 t). Color selec- 
tion methods such as Lyman- break (Steidel et al.lll995(l 
and narrow-band imaging techniques are providing statis- 
tical samples of very high redshift galaxies, allowing us to 
study large-scale structure at even earlier epochs. Lyman- 
break galaxies have correlation lengths as high as ro ~ 3 
Mpc even at z '--^ 3 — 4, and are thought to be as- 
sociated with (mildly) biased star-forming galaxies (e.g. 
Adelbergeil2000HOuchi et alJl200lHPorciani fc Giavahscol 
20021) . 



Studying clustering as a function of redshift and galaxy 
type may provide important constraints on some long- 
standing problems in cosmology concerning galaxy for- 
mation and evolution. For example, which of the galax- 
ies observed at high redshift are the progenitors of local 
galaxy populations, and which of the local galaxies host 
the remnant black holes that once powered high redshift 
active galactic nuclei (AGN)? Two common views on how 
structures observed at high redshifts may be related to 
structures observed today are represented by (i) the galaxy 
conse rvation model (e.g. iFrvl Il99(it iTegmark fc Peeblej 
Il998j) in which it is assumed that g a.laxies formed very 
early in a monolithic collapse (e.g. lEggen et alJ Il962l) 
and have evolved passively with a decreasing star for- 
mation rate si nce z ~ 2, and (ii) the hierarchical merg- 
ing model (e.g. IMo fc Whitelll99ql in which it is assumed 
that the most luminous galaxies formed more recently in 
massive dark matter haloes that have grown hierarchi- 
cally b y the merging of less massiv e galaxies and their 
haloes. iKauffmann fc Charlod l)l998() computed the evolu- 
tion of the observed K-band luminosity function for both 
the monolithic case and the hierarchical case, and found 
that by a redshift of ~ 1 these models differ greatly in the 
abundance of bright galaxies they predict. Likewise, the 
validity of these models can be tested by comparing pre- 
dictions for galaxy clustering from numerical simulations 

I II 1 

or ( semi-) anal vt hie t hcorv (e.g. Kauffmann et al. 1 999hl : 

iMoustakas fc Somerville 2002 : Mo fc White 2002. and ref- 
erences therein) with the observed clustering of a popu- 
lation of galaxies. In the case of pure monolithic collapse 
galaxy clustering is dictated by the evolution of galaxy 
bias under the rules of gravitational perturbation the- 
ory, but without the extra non-linear effects arising from 
galaxy mergers. Such a scenario can be thought of as a 



baseline model for the clustering of the matter as probed 
by galaxies situated in average mass haloes. However, 
in the hierarchical case the evolution of galaxy bias is 
much more complex, si nce galaxies are no longer conserved 
quantitities faauffma nn et all Il999bl). Comparing their 
observations to model predictions lDaddi et al ] ll2n0 find 
that such a scenario best explains the clustering evolution 
of massive ellipticals out to z ^ 1. 

Radio surveys can make an important contribution to 
this study: the use of magnitude-limited surveys for find- 
ing high redshift objects is usually a cumbersome task, 
while any flux density limited sample of radio sources con- 
tains objects at redshifts of z ~ — 5 ijPunlop fc Peacock 
E290). Powerful extra-galactic radio sources, or AGN in 
general, r esult from the fue lling of a supermassive black- 
hole (e.g. iReej 11984 Il99n|) . and there is evidence that 
the host galaxies of these high-redshift AGN are associ- 
ated with some of the most massive structures in the earlv 
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. Moreover, because powerful 



AGN were far more numerous at z ~ 1 — 2 than today, ra- 
dio surveys can be used to probe a population of massive 
galaxies in the epoch of galaxy formation. 

Despite initial concerns that any cosmological clus- 
tering of radio sources may be undetectable due to 
the relatively broad redshift distribution washing ou t 
the signal (e.g. 



Webster fc Pearsonlll977t lGriffithlll99l . 



iKooiman et al detected strong clusteri ng of bright 

radio sources in the 4.85 GHz 87GB survey. ICress et alJ 
(^1996.') made a thorough analysis of clustering at the 
mJv- level. Using the 1.4 GHz FIRST survey (see also 
iMagliocchetti et alJ Il998|) they obtained the first high- 
significance measurement of clustering from a deep radio 
sample, allowing them to investigate the separate con- 
tributions of b oth AGN and starburst galaxies (but see 
IWilman et al.l[2003. ). Further results on the statistics of 



radio sourc e clustering have been presented bv lLoan et alJ 
1I1997I) and lR,engelinlJ lll998V who based their analysis on 
the 4.85 GHz Parkes-MIT-NRAO survey and the 325 MHz 
WENNS survey, respectively. In high-resolution surveys 
such as FIRST, large radio sources can become resolved 
in several components, thereby spur iously contr i buting 
to t he cosmological clus t ering signal. ICress etal] lll996l) 
and ' Magliocchetti et alJ (|l998^ outlined the basic steps 
involved in separating the signal due to this effect from 
the true cosmological clustering, although the angular size 
distribution of radio sources at the mJy level is still largely 
unconstrained. 

Since the individual redshifts of the radio sources are 
generally not known, one usually only measures the two- 
dimensional clustering by means of the angular correlation 
function, w{9). However, the redshift distribution of the 
survey can be used to const rain rp. Using this so-ca lled 
Limber inversion technique llLimbed Il95a [Rub in, 1954t 
IPhillipps et al .ll978HPeeblesll98(t . radio sources from the 
above surveys are typically foui id to have rn ~ 5—15 
Mpc iRBngelink (,1998.) and iRengelink fc RottgeringI 
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l|1999|) pointed out that this broad range in rg measured 
can be explained by a scenario in which powerful radio 
sources have a larger vq than less powerful radio sources. 
This would be highly consistent with the mounting evi- 
dence that powerful radio galaxies are the high-redshift 
progenitors of local cD-galaxies residing in massive envi- 
ronments that are hence strongly clustered. Here, we will 
further explore the hypothesis of Rengelink et al. by inves- 
tigating the clustering of radio sources in a number of flux- 
limited subsamples t aken from the 1.4 GHz NRAO VLA 
Sky Survey (see also iBlake fc Wal]|l2002albh . the largest 
existing 1.4 GHz survey to date, containing ^ 1.8 x 10^ 
radio sources down to a flu x density limit of ~ 2.5 mJy 
at 45" (FWHM) resolution l)Condon et al.lll998(l . We also 
present new results on clustering using the latest release of 
the FIRST survey, carefully taking into account the con- 
tribution of multiple-component radio sources, which we 
found to be severely underestimated in earlier analyses. 

The outline of this article is as follows: in Sect. 2 we 
describe our methods for measuring the angular two-point 
correlation function. In Sect. 3 we describe the NVSS and 
FIRST radio surveys, and in Sect. 4 we present measure- 
ments of the angular clustering of the sources in these 
surveys and construct a simple model of the angular size 
distribution of radio sources. We derive an estimate of ro 
as a function of flux density limit in Sect. 5. In Sect. 6 we 
compare our results with the results found for other pop- 
ulations of galaxies taken from literature, and discuss how 
the combined measurements relate to current theories on 
galaxy formation and evolution. The main conclusions are 
summarized in Sect. 7. 



2. The angular correlation function 

The galaxy angular two-point correlation function, w(9), 
is defined as the excess probability, over that expected for 
a Poissonian distribution, of finding a gala xy at an angul ar 
distance 6 from a given other galaxy (e.g. |Peebiei[l983): 



foUow lReneelinkl l)l998|) and lWilman et all l)2003l) and use 

the Hamilton estimator 



SP = n[l + w{e)]6n, 



(1) 



where SP is the probablility, n is the mean surface den- 
sity and 6fl a surface area element. The angular two-point 
correlation function of a given sample of objects can be es- 
timated as follows. For each object, determine the angular 
distances to all other objects, then count the number of ob- 
jects in each angular distance interval, denoted by DD{6). 
As we want to calculate the excess probability of finding 
a galaxy at a certain distance from another galaxy due to 
clustering, we compare the observed distribution, DD{9), 
with the expected distribution of distances, RR{9), calcu- 
lated from large artificial catalogues of randomly placed 
sources. We note that several variants of w(6')-estimators 
exist in l i teratu re, of whic h the methods p r opose d by 
HamiltonI lll993l) and that of iLandv fc Szalavl l)l993() (see 



w{9) 



DD{9) ■ RR{9) 



{nD-l){nR-l)DR{9)-DR{9) 



- 1, 



(2) 



Blake fc Wal]ll2nn2b[ for application of this estimator to 



where no and ur are the number of sources in the data 
and random catalogues, respectively, and the numeri- 
cal factor 4n£)njj/(n£) — l)(ni^ — 1) normalizes the pair 
counts. This estimator additionally makes use of the cross- 
correlation between data and random catalogues, DR{9), 
to minimize effects due to large-scale fluctuations in the 
mean galaxy density. We estimate w{9) by averaging over 
the w{9) computed using 16 different random catalogues, 
each containing the same number of sources as the data 
catalogue to minimize the errors in DR{9) and RR{9) (a 
similar result can be obtained by constructing a single ran- 
dom catalogue that vastly exceeds the size of the data cat- 
alogue). Poissonian errors on the binned values of 'w{9) are 
estimated by 5w{9) = + w{9)]/ DD{9) . Alternatively, 
errors can be compu ted using the so-c alled bootstrap re- 
sampling method of iLing et al.l l)l986l) . In this method, 
the standard deviation in w{9) found among a large num- 
ber of pseudo-random resamples of the original dataset 
is used as a measure of the error in w{9). However, we 
found that fitting a model to w{9) (see Sect.0J using (i) 
Poissonian errors, and (ii) bootstrap errors gives results 
that are consistent within the errors of the fitted param- 
eters. Therefore, given the unprecedented volumes of the 
radio surveys we use the first method instead of the rela- 
tively expensive bootstrap technique. 

3. Survey descriptions and data selection 

3.1. The NRAO VLA Sky Survey 

The NRAO VLA Sky Survey (NVSS) is the largest ra- 
dio survey that currently exists at 1.4 GHz. It was con- 
structed between 1993 and 1998 JGondon et al.ll998^ ■ and 
covers - 10.3 sr of the sky north of J = -40° (~ 82% 
of the sky). Fig. ^ indicates the coverage of the NVSS. 
With a limiting fiux density of ~ 2.5 mJy {barms) and 
an angular resolution of 45" (FWHM), the NVSS con- 
tains about 1.8 X 10^ sources, and is consid ered to be 99% 
comp lete at a fiux density limit of 3.4 mJy l)Condon et alJ 
Il998|) . The NVSS is based on 217,446 snapshot observa- 
tions (of mostly 23 seconds) using the VLA in D- and 
DnC-configuration. These snapshots were then combined 
to produce a set of 4° x 4° datacubes containing Stokes 
I, Q, and U images. A source catalogue was extracted by 
fitting the images with multiple elliptical Gaussians. Since 
the angular resolution of the NVSS (9 w 45" FWHM) is 
well above the median angular size of extra-galactic radio 
sources {9 ~ 10 arcseconds), most sources in the catalogue 
are unresolved (> 95% for 3 < S'1.4 < 10 mJy). The main 
NVSS data products have been made publicly available 
for the use of the astronomical community, and can be 
obtained from the NRAO website^. 



NVSS) are generally considered to be the most robust. We ^ http://www.cv.nrao.edu/nvss/ 
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Fig. 1. AitofF map of the NVSS source density. Scales run from 2cr below (black) to 2(7 above the mean source density 
(white). The region of the galactic plane with |6| < 10° is indicated by solid lines. Besides the expected enhancement 
of the source density due to the large population of galactic radio sources, the NVSS catalogue suffers from large 
numbers of spurious sources around bright or extended sources (white regions), as well as an overall decrease in the 
source density below S = —10° (see the greyscale change at 6 = —10°). Sec text and Table^for details. 
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Fig. 2. The NVSS source density as a function of declination for various flux-limited sub-samples. Below ~ 10 mJy 
beam^^ the source density is non- uniform due to changes in the configuration of the VLA at (5 = —10° and 5 = +78° 
(dotted lines). 



3.2. NVSS data selection 

To optimize our catalogue for measuring the true cosmo- 
logical clustering of radio sources, we have carried out 
a detailed examination of the NVSS source catalogue 
to identify and correct regions that may spuriously 
contribute to w{9) : 

(i) The edge of the survey just a few arcminutes south 
of S = —40° follows an irregular pattern with right 
ascension. We select the region 6 > —40° to ensure that 
the boundary of the survey is straight. 

(ii) The survey area is known to contain six hexagonal 
gaps due to missing snapshot observations that we 
masked from the catalogue by excluding rectangular 
regions of 2° x 2° fully covering each gap. The regions are 
listed in Table □ 

(iii) We constructed a map of the NVSS source density as 
a function of position on the sky by applying an equal-area 



projection to the catalogue and plotting filled contours 
of the number of objects in 1° x 1° non-overlapping cells 
covering the survey area. This map is shown in Fig. ^ 
The scaling of the greyscale was chosen so that under- 
dense regions of 2a below the mean density are black, 
and overdense regions of 2a above the mean are white. 
Radio emission from the region of the galactic plane, as 
evidenced by a continuous chain of large white areas in 
Fig. n is dominated by the large population of galactic 
radio sources that consists mostly of supernova remnants 
and HII regions. In Fig. 0] we plot the rms- noise level as a 
function of galactic latitude, where the rms-noise level in 
each latitude bin is the average of the locally determined 
rms-noise values listed for every source entry in the NVSS 
catalogue. The rms-noise level is found to peak at 6 = 0° 
due to the overcrowding of galactic sources, but falls off 
to a relatively constant level of ^ 0.48 mJy beam"-'^ for 
|5| > 10°. We decided to exclude the region of the galactic 
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Fig. 4. The rms-noise level as a function of galactic lati- 
tude. The average rms-noise level of the survey is ^ 0.48 
mJy beam~^. Dotted lines enclose the region \b\ < 10°. 



plane that is bounded by |5| — 10°, which was chosen 
so that the large overdense regions in Fig. ^ are all fully 
masked and the rms-noise is at a relatively constant level. 

(iv) Further inspection of Fig. ^ reveals that some regions 
are associated with a significant increase in the local 
source density. From contour maps of these areas it was 
found that bright and/or extended sources are sometimes 
accompanied by significant numbers of spurious sources 
due to a side-effect of the fitting algorithm used to extract 
the sources, and, in some cases, due to side-lobe con- 
tamination. From the catalogue we excluded rectangular 
regions of mostly 1° x 1° in size centered on each of these 
sources (larger regions of up to 2° x 2° were required in 
some cases). The excluded regions are listed in Tabled 
No regions of > 2a underdensities were found. 

(v) Most of the NVSS observations were conducted 
using the VLA in D-configuration, but the regions 
S < —10° and S > -f78° were observed using the hybrid 
DnC-configuration to counterbalance projection effects 




S > 3 mJy 

S > 5 mJy 
S>7mJy 
S > 10 mJy 

S > 20 mJy 



10 20 30 40 
Declination (°) 



50 60 



Fig. 5. The FIRST source density as a function of decli- 
nation for various limiting flux densities. 



which result from foreshortening of the north-south 
uw-coverage range. Fig. |21 shows the NVSS source density 
as a function of declination for various flux-limited 
sub-samples. Below the flux density limit of 10 mJy, the 
use of the DnC-conflguration has caused a signiflcant 
decrease in sensitivity leading to a drop in the source 
density of > 10 percent (see also Fig. QJ. As this will 
inevitably cause spurious signal in the angular two-point 
correlation function, we selected only the regions observed 
in D-conflguration for measuring w{6) below flux density 
limits of 10 mJy. 

Table|211ists the final regions and the number of sources 
in them for various fiux density limited subsamples. 



3.3. The FIRST Survey 

The FIRST (Fain t Images of the Radi o Sky at Twenty cen- 
timeters) survey l|Becker et al.lll995|) is another 1.4 GHz 
VLA survey, which was started in 1993 and is still under 
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Table 1. Regions of the NVSS catalogue that were masked because of missing snapshot observations and overdense 
regions associated with bright or extended sources. Overdense regions at |6| < 10° are not listed here since we excluded 
this area from the catalogue as a whole. 



RA (J2000) DEC (J2000) Remark 



15''38'"00^ 


- 15'^ 43^00" 


-05°00'00" 


- -06°00'00" 


Missing snapshot 


09'' 54'" 00^ 


- 10'' 00"" 00" 


-11°45'00" 


- -12°45'00" 


Missing snapshot 


09'' 54'" 00^ 


- 10'' 00'" 00" 


-25°00'00" 


- -26°00'00" 


Missing snapshot 


04'' 25'" 00" 


- 04'' 30"" 00" 


-37°45'00" 


- -38°45'00" 


Missing snapshot 


18''17"'00" 


- 18''22'"00" 


-16°00'00" 


- -17°00'00" 


Missing snapshot 


18''02"'00" 


- 18''07'"00" 


-23°30'00" 


- -24°30'00" 


Missing snapshot 


0l''34"'00" 


- 0l''42'"00" 


+32°30'00" 


- +33°50'00" 


3C 48 




- 03''21'"48" 


+40°00'42" 


- +43°00'42" 


Perseus A 


03'' 17"' 00" 


- 03'' 27"" 00" 


-36°20'00" 


- -38°20'00" 


Fornax A 


04'' 35™ 05" 


- 04'' 39'" 05" 


+29°10'12" 


- +30° 10' 12" 


3C 123 


05'' 18^00" 


- 05''26"'00" 


-35°40'00" 


- -37°20'00" 


PKS 0521-36 


05''31"'17" 


- 05''39'"17" 


-06°23'00" 


- -04°23'00" 


M42 


05''38"'00" 


- 05''46'"00" 


-01°20'00" 


- -02°40'00" 


3C 147.1 


05''40"'36" 


- 05'' 44™ 36" 


+49°21'07" 


- +50°21'07" 


3C 147 


05''52"'00" 


- 05''56™00" 


-04°20'00" 


- -05°40'00" 


TXS 0549-051 


07''05"'00" 


- 07''25™00" 


+74°20'00" 


- +75°20'00" 


3C 173.1 


09''15"'05" 


- 09''21™05" 


-12°50'24" 


- -11°20'24" 


Hydra A 


12''16"'00" 


- 12''23™00" 


+05°00'00" 


~ +06°30'00" 


NGC 4261 


12''26"'07" 


- 12''32™07" 


+01°18'00" 


- +02°48'00" 


3C 273 


12'' 26™ 50" 


- 12''34™50" 


+11°23'24" 


- +13°23'24" 


M87 


13''11"'00" 


- 13''15™00" 


-22°30'00" 


- -21°30'00" 


MRC 1309-216 


13''21'"00" 


- 13'' 27™ 00" 


+31°20'00" 


- +32°20'00" 


NGC 5127 


14''07"'00" 


- 14''15™00" 


+51°40'00" 


- +52°40'00" 


3C 295 


16''49"'11" 


- 16''53™11" 


+04° 29' 24" 


- +05°29'24" 


Hercules A 


17'' 18^00" 


- 17''24™00" 


-01°40'00" 


- -00°20'00" 


3C 353 


18''26"'00" 


- 18''34™00" 


+48°00'00" 


- +49°40'00" 


3C 380 


19'' 22^00" 


- 19''26™00" 


-28°45'00" 


- -29°45'00" 


TXS 1921-293 



construction. Using the VLA in B-configuration it will 
ultimately cover ~ 10, 000 square degrees of the north- 
ern Galactic cap, matching the survey area of the Sloan 
Digital Sky Survey. Given the large coverage of FIRST, 
its sensitivity is unprecedented: with a limiting flux den- 
sity of ^ 1 mjy (barms) and an angular resolution of 5'.'4 
(FWHM) the catalogue contains about 100 sources per 
squar e degree with a co mpleteness level of ~ 95% at 2 
mJy l)Becker et al.lll995|) . 

We have obtained the publicly available 2001 October 
15 version of the source catalogue^, which has been de- 
rived from the 1993 through 2001 observations, and cov- 
ers about 8,565 square degrees of the sky. About 4% of 
the 771,076 sources in the catalogue are flagged as possi- 
ble side-lobes, which we exclude from the catalogue. We 
set the lower flux density limit of the catalogue to 3 mJy, 
the limiting flux density of the NVSS survey. Finally, we 
select the regions +2° < ^ < +20° and 9'' < a < 16'', 
+20° < ^ < +55° and 8'' < a < 17'' from the catalogue, 
by requiring a relatively uniform source density and a sim- 
ple geometric form. This area covers ~ 5, 538 square de- 
grees and contains 188,885 sources. As for the NVSS, we 
construct a map of the FIRST surface density (Fig.OJ. and 

^ http://sundog.stsci.edu/ 



plot the source density as a function of declination (Fig. 
ISJ . For the selected region we found no suspicious features 
in the catalogue. The number of sources in various FIRST 
subsamples are listed in Table El 

4. The angular clustering of radio sources 

4.1. The angular correlation function of S > 10 mJy 
NVSS sources 

Following the procedures described in Sect.|21we compute 
■w{9) for the 5 > 10 mJy NVSS subsample. Distances be- 
tween data and/or random positions are initially measured 
in bins of 0f5, and rebinncd in bins of constant logarith- 
mic spacing to analyse the data. We fit the data using a 
weighted ^^-minimization routine, and we determine the 
la errors from the covariance matrix. 

The results are shown in Fig. We find that two 
power-laws are needed to describe the full range of our 
measurements. Fitting the data with a p ower-law angu lar 
correlation function w{9) = AO'^-'f (e.g. IPeebleslllisol) at 
angular scales of < 6' gives a slope of 7 = 4.4 ± 0.2, 
while at 6* > 6' we find a slope of 7 = 1.7 ± 0.1. The 
latter value is consistent with the slope of the empirical 
power-law of 7 ~ 1.8 found for the cosmological cluster- 
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NVSS 






FIRST 


Slow Region 


Sources 


Slow 


Region Sources 


3 mjy 


+ 10° < & < +45°, -5° < <5 < +70° 


210,530 


3 mJy 


+2° <5 < +20° and 9" < a < 16" 188,885 
+20° <S< +55° and 8'' < a < 17*" 


5 mJy 


b 


> 10°, -10° < <5 < +78° 


507,608 


5 mJy 


" " 124,974 


7 mJy 


b 


> 10°, -5° < 5 < +70° 


351,079 


7 mJy 


" " 94,099 


10 mJy 


b 


> 10° 


433,951 


10 mJy 


" " 68,560 


20 mJy 


b 


> 10° 


242,599 






30 mJy 


b 


> 10° 


165,459 






40 mJy 


b 


> 10° 


123,769 






50 mJy 


b 


> 10° 


97,753 






60 mJy 


b 


> 10° 


79,738 






80 mJy 


b 


> 10° 


56,903 






100 mJy 


b 


> 10° 


43,294 






200 mJy 


b 


> 10° 


17,015 








0.0001 
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Fig. 6. The angular two-point correlation function of 5 > Fig. 7. The angular two-point correlation function of 5 > 
10 mJy NVSS sources. The power-law fits described in the 3 mJy FIRST sources. The power-law fits described in the 



text are indicated. 



text are indicated. 



ing of objects ranging from normal galaxies to clusters 
Csee iBahcall fc SoneiralHoS^. for a review). However, at 
small angular scales the power-law is much steeper, pre- 
sumably caused by the enhancement of DD{9) due to 
the decomposition of large radio galaxies into their sep- 
arate radio components (see Sect. l4.2l fc Sect. 14.41 see also 
Blake & Wall 20G2a). If we fit the data simultaneously 
with a double power-law correlation function of the form 
w{e) = B9^"'^^ + A6'i-^^ with fixed slopes of 75 = 4.4 
and 7A = 1.8, we find amplitudes of B = (1.5±0.2) x 10"^ 
and A = (1.0 ± 0.2) x lO'^. The double power-law fit is 
indicated in Fig. El 

4.2. The effect of multiple component radio sources 
and the angular correlation function of FIRST 

Alth ough the median ang ular size of radio sources is ^ 10" 
(e.g. ICondon et ahlHoOsli . radio sources can have sizes of 
up to several arcminutes. At angular scales comparable 



to the size of these large radio galaxies, the true cosmo- 
logical w{9) can become confused or even dominated by 
resolving these galaxies into their various radio compo- 
nents, such as lobes, hot spots and cores. The angular scale 
at which the size distribution of radio galaxies begins to 
dominate wlO) is indicated by the clear break around 6'. 
Earlier studies attempted to correct w{9) for the contribu- 
tion of multi-component radio sources by me ca,ns of com- 
ponent combining algorithms. For example, ICress et alJ 
( 1996) calculated the angular correlation function for the 
FIRST survey considering all sources within 1.'2 of each 
other as a single source. The analys i s of th e FIRST data 
was repeated bv lMagliocchetti et alJ l)l998l) . who removed 
double sou rces using a i i algo rithm based on the 9 oc VS 
relation of lOort et al ] lll987l) and flux ratio statistics of 
the components of genuine doubles. They found values of 
7 = 2.5 ± 0.1, and A = (1.0 ± 0.1) x IQ-^ for fiux den- 
sity limits between 3 and 10 mJy. Comparing their results 
to our measurement for the NVSS presented in Fig. El we 
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Fig. 8. Angular correlation functions for the flux density 
intervals 10 < S* < 40 mJy and S > 200 mJy. The power- 
law fits to the data described in the text are overplotted. 
Because of an unexplained 'bump' in the S > 200 niJy 
signal at 0.1 < 6* < 0.3 (connected points), the small- and 
large-scale correlation functions were fitted separately over 
the ranges 6 < 0.1 and 9 > 0.3, respectively. 



conclude that despite the efforts of these authors it is likely 
that a residual contribution from large radio galaxies re- 
mained. Fitting the data over the whole range of 9 with 
a single power-law explains the apparently high value of 
7 ~ 2.5 reported for the clustering of FIRST radio sources. 

Here, we present new measurements from the FIRST 
survey . Our reasons for repeating the work of ICress et alJ 
l) 19961) and Maglio cchetti et al. (199"^ are threefold. 
Firstly, the FIRST catalogue has almost doubled in size, 
enabling a better statistical measure of 'w{9). Secondly, the 
clear break found in the angular correlation function of the 
NVSS enabled us to isolate the signal due to true cluster- 
ing from the signal due to the size distribution of radio 
galaxies. A similar analysis can be applied to the FIRST 
data. Thirdly, we found large-scale gradients in the NVSS 
source density below a flux density limit of 10 mJy (see 
Sect. E21l- The FIRST data can be used to verify and 
complement the results from the NVSS for 3-10 mJy. 

In Fig. 13 we present our measurements for the angular 
correlation function from the S* > 3 mJy FIRST subsam- 
ple. As for the NVSS, we see a clear break in 'w{9) due 
to the presence of multi-component radio sources. Fitting 
the measurements with our double power-law model yields 
7s = 4.1 ± 0.2 and B = (2.7 ± 0.3) x 10^^, and 7^ = 
1.9±0.2 and A = (1.0±0.3) x 10~3. Note that the break in 
w{9) in this sample occurs at 9 ^ A' compared to 9 ^ 6' for 
5* > 10 mJy in NVSS (see Fig-. IHll. iBla.ke fc Walll l|2002afl 
show that this is due to a 1/cr dependency [a being the sur- 
face density of radio sources) of the amplitude of ■w{9) at 
small angular scales, simply because the weight of pair- 
counts due to large radio galaxies increases as the surface 
density decreases (see their equation 4). 



1 0.010 

I 
< 



0.001 




10 100 
1.4 GHz flux density limit (mJy) 

Fig. 9. The amplitude of the cosmological angular correla- 
tion function (7 = 1.8) of NVSS and FIRST as a function 
of 1.4 GHz flux density limit. For comparison, we have in- 
dicated the results for the WENSS and GB6 surveys from 
Rengelink, (,19.9E.) and Rengelink & R.ottgerina (1999). 

We conclude that the cosmological 'w{9) of 5 > 10 
mJy NVSS sources and S > 3 mJy FIRST sources, as 
determined by our analysis, are consistent with having 
the canonical clustering power-law slope of 7 ~ 1.8, and 
an amplitude of A ~ 1 x lO^'^. 

4.3. w{9) as a function of flux density limit 

To investigate angular clustering as a function of flux den- 
sity hmit, we calculate w{9) for aU NVSS and FIRST sub- 
samples listed in Table 12 We obtain the amplitudes of 
'w{9) by fitting the data with the double power-law model 
w{9) = B9^-'<'' + A9^-''^, fixing the slopes at 75 = 4.4 
and = 1-8. However, because the signal for the S > 200 
mJy subsample is affected by a 'bump' at 6* 0?2 (see 
Fig. IHl) , we obtained the amplitudes for this subsample by 
fitting the small- and large-scale correlation functions sep- 
arately with power-laws w{9) = B9~^'^ for 6* < 0?1 and 
w{9) — A9^^-^ for 9 > 0?3, respectively. The measured 
amplitudes and their la errors are listed in Table 01 The 
values of both B and A are found to increase with increas- 
ing fiux density limit of the subsamples. The increase in B 
can be explained by the 1 / cr-dependency of the small-scale 
correlation function that is dominated by double or mul- 
tiple component radio sources (see Sect. 14.2(1 . From this 
point onward, we will be only concerned with the ampli- 
tude A that is believed to be dominated by the true cosmo- 
logical clustering. In Fig.|51we have plotted the amplitude 
of the cosmological 1^(9) as a function of fiux density limit. 
For comparison, we have indicated the result s from the 
325 MHz WENSS and 4850 MHz GB 6 surveys l|Rengehnkl 
Il998t iRengelink fc Rottgerindll999() by extrapolating to 
1.4 GHz using a power law spectrum. Si, oc i'^", with 
spectral index a ~ 0.8. Between 3 and 40 mJy the ampli- 
tude is approximately constant within the errors and has 
an (unweighted) average of ~ 1.2 x 10~^. From 50 — 100 
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Table 3. Amplitudes and la errors of the double power- 
law correlation function w{6) = BO^'^ '^ + A9^-^ as a func- 
tion of flux density limit. 





NVSS 


FIRST 




10" X B 


10"* X A 


10'' X B 


10^ X A 


3 mjy 


0.7 ±0.1 


0.8 ±0.1 


0.7 ±0.1 


1.2 ±0.1 


5 mJy 


1.0 ±0.1 


1.0 ±0.1 


1.0 ±0.1 


1.4 ±0.2 


7 mJy 


1.2 ±0.1 


0.9 ±0.2 


1.2 ±0.1 


1.3 ±0.3 


10 mJy 


1.5 ±0.1 


1.0 ±0.2 


1.4±0.1 


1.9 ±0.3 


20 mJy 


2.6 ±0.2 


0.9 ±0.2 






30 mJy 


4.1 ± 0.1 


1.0 ± 0.2 






40 mJy 


4.4 ±0.1 


1.3 ±0.5 






50 mJy 


8.3 ±0.2 


1.9 ±0.8 






60 mJy 


6.8 ± 3.0 


1.7 ±0.4 






80 mJy 


8.4 ±0.2 


1.9 ±0.7 






100 mJy 


19 ± 1 


2.2 ± 0.7 






200 mJy 


30 ± 1 


6.6 ± 1.8 







mJy the amplitude is ~ 2x higher, and it has increased by 
another factor of ~ 2 — 3 at 200 mJy. These measurements 
indicate a trend of increasing clustering amplitude with 
increasing flux density limit. However, one has to keep in 
mind that the sources in the brighter subsamples are also 
included in the subsamples with lower limiting flux den- 
sities. Therefore, we also compute w{6) for sources that 
lie in the flux interval 10 < S" < 40 mJy. The results are 
shown in Fig.|Sltogether with found for S > 200 mJy. 
The amplitude A — (6.6± 1.8) x lO^'^ that we measure for 
S > 200 mJy is significantly higher than the amplitudes 
measured at low er fl ux densities. This is cons istent with 
iReneehnkl l)l998() and lReneelink fc Rotteerind figQa) who 
found A = (11.5 ± 3.5) x IQ-^ for Sia > 160 mJy in the 
GB6 survey and 'Loa n et alJ l)l997l) who estimated that A 
has a value between 0.005 and 0.015 for S'1.4 > 100 - 270 
mJy from the combined 87GB and PMN surveys (Fig. ^ . 
W e would like to ma ke th e following remarks: 

(i) iRengelinkI l)l998|) and iRengeUnk fc Rotteerin j l)l999(l 
measured w{9) from WENSS and GB6 by excluding the 
first 5' and 10', respectively. We have used our routines 
to measure w{6) for their catalogues as well (not shown 
here). The amplitudes and slopes we find are consistent 
with their values, and we find no evidence for a contribu- 
tion of multi-component sources at the smallest angular 
scales allowed by these surveys. 

(ii) Below 10 mJy the amplitudes for the NVSS and 
FIRST data are consistent with A ~ 1.1 x 10~^. However, 
at 10 mJy the amplitude is '--^ 2x higher for FIRST than 
for the NVSS. This is curious since the NVSS and FIRST 
surveys probe radio so urces at exactly the same frequency. 
iBlake fc Wall ll2002bl) give a very nice demonstration (see 
their Fig. 3) of the most probable cause. The resolution of 
FIRST is ten times higher than that of NVSS, and there- 
fore the average flux density of a single NVSS source is 
only equal to the sum of all its possibly resolved compo- 
nents in FIRST. Sources that appear in NVSS with in- 



tegrated fluxes just above a given flux density limit can 
thus be missed in FIRST. Therefore, we consider NVSS 
to be more optimal than FIRST for measuring the clus- 
tering of extra-galactic radio sources. Furthermore, if we 
compute w{9) for only those NVSS sources that lie in 
the region covered by FIRST, we flnd an amplitude of 
A = (1.7 ± 0.3) X 10"''. This is consistent with the re- 
sults found for the 10 mJy FIRST sample, suggesting that 
cosmic variance of clustering may be an additional factor 
contributing to the difference in amplitudes measured for 
the total NVSS area and FIRST. Future work might show 
that the region covered by FIRST is especially rich in 
large-scale structures. 

(iii) In the 200 mJy subsample we find an unexpected in- 
crease in the correlation signal at 6* « 0?2 (indicated by 
the connected points i n Fig. Et- W e inves tigate two pos- 
sibilities. (1) Sidelobes: ICress et all l|l996|) found a bump 
in w(0?l) for S" > 3 mJy sources in FIRST, and found 
that it was caused by sidelobe contamination. However, 
if sidelobes are responsible for boosting the correlation 
function at 6* ~ 0? 2 in the S > 200 mJy NVSS sample, 
these sidelobes themselves also must have minimum peak 
fluxes of 200 mJy. It is highly unlikely that such bright 
sidelobes have found their way into the NVSS catalogue, 
without being masked in Sect. 13.21 Also, we have visually 
inspected the contour maps of several tens of source pairs 
{S > 500 mJy) that contribute to w{0) at 9 ^ 0?2. In all 
cases the pairs consisted of unresolved peaks without signs 
of diffuse, extended emission or side-lobe contamination. 
(2) Radio galaxies with large angular sizes: the position 
of the bump near the break in w{9) suggests that it may 
somehow be relat ed to the size distribution. Conveniently, 
iLara et al] l)200l[) have constructed a sample of 84 large 
angular size {9 > 4') radio galaxies from the NVSS at 
S > ±60° and a total integrated flux density of > 100 
mJy. Candidates were pre-selected by visual inspection of 
the NVSS maps, and conflrmed or rejected following ob- 
servations at higher resolution. If the bump is caused by 
~ 12'-sized radio galaxies, then giv en the 2-Mpc linear 
size c utoff of large radio galaxies fsee lSchoenmakers et alJ 
I2OOII) . these galaxies must lie at z < 0.1. It is unlikely 
that such a large, relatively nearby source with, among 
other emission, two radio components each with a peak 
flux of > 200 mJy would have been missed by their selec- 
tion criteria. iLara et all l)200lt) determined angular sizes 
by either measuring the maximum distance between 3a 
contours, or by the distance between peaks at the source 
extremes. Also, sizes were measured along the 'spine' of 
a source if significant curvature was present. To investi- 
gate how many of these sources could actually contribute 
to w{9) at ^ 12' we r edetermine the angular sizes of the 
sources of iLara et al.l |2001). We find that none of these 
sources consists of > 2 components of > 200 mJy of ^ 12' 
separation. On the other hand, if we extrapolate the clus- 
tering power-law derived at larger scales to 6* = 0? 2 we 
find that the bump translates into ~ lOx the number of 
pairs expected. Even allowing for the much larger area of 
NVSS, the possibility that the bump is caused by large 
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radio galaxies as in the sample of iLara et alJ 1)200 ]|) is 
therefore unlikely. 

Unfortunately, the exact origin of this feature remains 
unclear. We realize, however, that this bump is situated at 
a crucial angular scale for our measurements. Therefore, 
we have obtained the amplitudes B and A by fitting 
w{9) on both sides of the bump with a single power-law. 
Under the condition that the effect that causes the bump 
is not responsible for enhancing w{8) at > 0?3, this will 
enable us to derive an estimate for the amplitude for the 
cosmological clustering. At ^ > 0?3 •w{0) is consistent with 
the classical 7 = 1.8 power-law clustering model. 

4.4. Modelling the angular size distribution of radio 
galaxies 

4.4.1. The model 

The steepening of the slope oiw{9) at small angular scales 
is presumably related to multi-component sources spuri- 
ously enhancing the true clustering pair counts at small 
9. To demonstrate the reality of this assumption, we cre- 
ate a simple model for the angular size distribution of ra- 
dio galaxies in the NVSS, that is able to account for this 
extra signal contributing to w{9). We model the physi- 
cal size distribution of sources in our S* > 10 mJy NVSS 
sample, and use their redshift distribution to obtain the 
angular size distribution. Because we know the angular 
resolution of the NVSS, this model can then be used to 
estimate the fraction of sources likely to be resolved. It is 
essential to separate sources that are resolved into a sin- 
gle, elongated object from sources that are resolved into 
a number of components, since only the latter would pro- 
duce extra pair counts. Here, we assume that the majority 
of surplus pair counts arise from resolving the two edge- 
brightened radio lobe s of FRII-type radio galaxies (see 
[Panaroff & Riley Il974l) , and we es timate that the fraction 
of FRIIs at 10 mJy is - 40% from lWall & JacksonI l|l997t) 
(assuming a spectral index of a = 0.8 to extrapolate to 
1.4 GHz). 

Several groups have investigated the median physi- 
cal sizes of FRII radio galaxies as a function of redshift 
and radio luminosity by parameterizing the linear size as 
-D oc (1 -f z) ~'^P™, where P is t he radio luminosity (for 
a review see iBlundel l et al]ll999|) . Results using different 
sam ples of rad io galaxies vary from no size evolution at all 
Ce.g. lNilssonet al.lll993ll. to size evolu tion depending only 
on redshift (e.e. ^Kaoahi et alJ ll987t) . and size evolution 
depending on both redshift an d luminosity with contradic- 
torv results (e.g. lOort et a,lJll'98i iBarthe! fc Mil^Jigsi 
ISinga]|ll99j) . We use the results of iNeeser et all l)l995|) 
who found the following linear size-redshift relation from 
a spectroscopically complete sample of FRII radio galax- 
ies: 

L»cx (l + z)^^-^^"-^ (for 17m = 1 and = 0), (3) 

and remark that no intrinsic correlation was found be- 
tween D and P (P™ ~ 1 with m = 0.06 ± 0.09). This 



observed linear-size evolution may be related to evolution 
of the confining intergalactic medium, or to evolution of 
the radio galaxy itself, but the exact und erlying physical 
mechanism is unknown fsee lNeeser et al.lil995.V 

For the purpose of our model, we place simulated 
sources in small redshift intervals (Az = 0.01) in the 
range < z < 5, and assume that their mean phys- 
ical size evolves with redshift according to equation 
We set the total number of input sources equal to the 
estimated number of S" > 10 mJy FRIIs in our NVSS 
sample (~ 40%o/434, 000), and calculate the number of 
sources in each redshift interval from the redshift distri- 
bution, N{z), using the formalism of IPunlop fc Peacockl 
(1990) (sec Sect. El for details). We then assume that in 
each redshift interval sizes are normally distributed. We 
take a mean size of 500 kpc and a standard deviation 
of 250 kpc at z ~ 0, chosen so that the resulting phys- 
ical size distribution roughly resembles the distribution 
of projected linear s izes versus redshift as it is given by 
iBlundell et alJ ljl999l) for three complete samples of FRII 
radio galaxies from the 3C, 6C, and 7C radio surveys. 
The resulting physical size distribution is shown in Fig. 
[TUl where we plot filled contours of the source density in 
the linear size-redshift plane to illustrate the underlying 
redshift distribution. We have also indicated the minimum 
physical size that is theoretically required for a source to 
become resolved as a function of redshift, given by the 
NVSS resolution of 45" (FWHM). We would hke to re- 
mark at this point that the distribution of sizes in our 
model beyond redshifts of 2 ~ 3 should not be taken too 
seriously as it is based on a straight extrapolation from 
measurements made at redshifts < z < 2, and does 
not take into account the fact that at these high redshifts 
most sources will be extremely young and are thus likely 
to be very small. However, as can be seen from Fig. 1101 our 
modeled size distribution falls below the NVSS resolution 
already at z ~ 1. Taking smaller sizes at higher redshifts 
will have no effect on the modeled size distribution of re- 
solved sources that we want to derive here. 

Assuming VIm = 1 we calculate the angular size distri- 
bution associated with our model. We construct 10 such 
models, and average them to get our final model of the 
angular size distribution of the sample. This model is pre- 
sented in Fig.^] Although the m e an an gular size is ^^10" 
in agreement with ICondon et alJ ((1998') , sizes are found 
to extend up to several arcminutes beyond the resolution 
of the NVSS (indicated by the dotted line). 

4.4.2. Results 

We now compare the number of surplus pairs expected 
from resolved FRII sources in the model, DDmod{9), to 
the actually measured pair counts at angular scales oi9 < 
6'. At these scales, the measured pair counts consist of 
both pair counts due to clustering and pair counts due to 
doubles, so 

DDtot{9)=DDgai{9)+ DDdu (9). (4) 
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Fig. 10. The modeled physical size distribution of S" > 10 
mJy FRII radio galaxies in the NVSS catalogue. The 
source density in the linear size-redshift plane is indicated 
by contours to illustrate the underlying redshift distribu- 
tion (darker greyscales indicate higher densities). Sources 
lying above the line can, in principle, be resolved given 
the angular resolution of the NVSS of 45"(FWHM). 

To extract DDdbi{d) from the total counts, DDtot{d), we 
calculate DDgai{0) by assuming that the galaxy angu- 
lar correlation function as measured above the break in 
i«(6') can be extrapolated to angular scales of < 6': 

Wgai{e) = 1.0 X lO-30-° » = DDgaiiO) ■ F{e) - 1, (5) 

where F{e) = 4 • RR{9)/[DR{e)]^ , the part of the 
Hamilton estimator that is relatively independent of the 
presence of doubles. Since we now know both DDtot{0) 
and DDgai{0), we can subtract them to get a measure of 
the counts arising from the double sources: DDdbi{0). The 
final step is to rebin the modeled number of pair separa- 
tions DD„iod{0) in order to match the binning scheme of 
DDdbi{d)- Fig. ll2lshows the ratio of the observed doubles 
to the modeled doubles per distance interval. The errors 
in the observed counts are estimated from the la-error 
in the amplitude of 'w{9). The errors in the modeled pair 
counts are estimated by allowing a 10% error in the esti- 
mated fraction of FRIIs in the NVSS. We conclude that: 
a model in which the small-seale angular correlation func- 
tion steepens due to resolving FRII radio galaxies into two 
distinct knots of radio emission is in good agreement with 
the measurements presented in Fig. O 

Several remarks that can be made are the following: 
(i) The size distribution of radio sources at the mJy level is 
still la rgely unconstrained. Recently, however. iLara et al.l 
l)200ll) presented a new sample of large radio galaxies 
(LRGs) selected from the NVSS. In the region 5 > +60° 
they found ~80 radio galaxies with apparent angular sizes 
larger than 4' and total flux density greater than 100 mJy. 
If we roughly extrapolate our model to their sensitivity 
and correct for the area we successfully predict the num- 
ber of FRIIs in the range 4' < 6* < 6^ However, in this 
interval one third of the sample of iLara et Til ll200lh con- 
sists of FRIs, while the model only uses FRIIs to estimate 
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Fig. 11. The angular size distribution for FRII radio 
galaxies in the NVSS calculated from the modeled phys- 
ical size distribution (assuming CIm = 1). The number of 
input sources was chosen to match the predicted number 
of FRIIs in the 5 > 10 mJy subsample. The binsize is 1". 

the number of surplus pairs expected. The model could 
be refined by decreasing the fraction of resolved FRIIs to 
also allow a contribution from large FRIs. 

(ii) The model allows objects to be either single or dou- 
ble sources, although visual inspection of NVSS contour 
maps shows that sources are sometimes split into three 
or even more components. Therefore, we may expect an 
extra amount of spurious pair counts on top of the counts 
due to classical double radio sources. This may become 
increasingly important with increasing flux density limit. 

(iii) The model predicts a fraction of resolved sources in 
NVSS of ^0.07, i n rough agreement w ith the value of 
-0.05 predicted bv lCondon et alJ l)l998|) . 

The simple model allows us to explore the general 
relations between the physical size distribution of ra- 
dio galaxies and w{9) at small angular scales. Although 
our crude method is successful in reproducing the obser- 
vations, it relies on a number of assumptions that are 
not easily verified from the data currently in literature. 
Radio sources come in a wide variety of sizes ranging 
from < 1 kpc for the class of gigahertz peaked spec- 
trum sources (GPS), to 1-20 kpc for the compact steep 
spectrum sources (CSS), > 20 kpc for FRI- and FRII- 
type radio g alaxie s, and > 1 Mpc fo r giant radio ga laxies 
llFanti et a1.lll flflOHO'Dea et a.ljFflmllPilundell et a.lJlTflfl^ 
ISchoenmakers et al.ll200ll) . Evidently, the distribution of 
linear sizes of radio sources are very complex, and will re- 
main an important subject for future studies. As we have 
shown, the angular correlation function can be used to put 
constraints on the size distribution of large radio galaxies. 
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9 (arcmin) 

Fig. 12. The ratio of observed doubles to modeled doubles 
per distance interval. The angular resolution of the NVSS 
is indicated by the dotted line. 



where V(z) is the comoving volume, $i(P, z) is the 
model RLF, x{z) the comoving distance, S the limiting 
flux density of the subsample, and v the frequency of 
FIRST/NVSS. We note that N{z) is independent of cos- 
mology as long as the calculations are carried out in the 
cosmology used to construct the RTFs (i.e. Hm = 1-0 and 
Ho = 50 km s-i Mpc"^). 

Figs . 1 1 31 and 1 1 41 show the rcdshift distributions for S > 
10 mJy and S > 100 mJy, respectively. We calculate the 
average of the six different models (indicated by the solid 
curve), which will be our best estimate of N{z) use in 
the analysis below (the same method was used for the 
N{z) applied to the model of the angular size distribution 
described in Sect. l4.4|l . It is important to keep in mind that 
the functional form of N(z) remains virtually unchanged 
from 3 — 200 mJy. Over this range in flux densities the 
RLFs represent a broad redshift distribution with a peak 
around z ~ 1, indicating the very large median redshift 
that is generally probed by radio surveys. 



However, perhaps more ideal would be to make a statis- 
tical redshift sample of all radio source pairs within some 
angular distance interval, and then take high resolution 
radio observations to constrain the numbers of intrinsic 
doubles in that sample. 

5. The spatial clustering of NVSS souces 

5.1. The redshift distribution 

At the mJy level and higher it is stand ard practice to com- 
pute r edshift distributions using the 'Punlop fc Peacocl] 
IjlQQOl) radio luminosity functions (RLFs). These authors 
have constructed a range of model luminosity functions us- 
ing spectroscopically complete samples from several radio 
surveys at different frequencies. Using a free-form mod- 
elling approach they found a number of smooth func- 
tions that were consistent with the data. In addition, 
they attempted two models of a more physical nature by 
assuming pure luminosity evolution (PLE) and luminos- 
ity/density evolution (LDE) to describe the RLF. The to- 
tal ensemble is expected to agree well at those luminosities 
and frequencies at which they are best constrained by the 
data, while uncertainties in the extrapolation of each of 
these models to those regions that are less constrained 
by the data may be reduced by taking the ensemble as a 
whole. We compute redshift distributions, N{z), for each 
flux-limited subsample using the free-form models 1 — 4 
and the PLE/LDE models for the combined population of 
flat (a = 0, S,y (X z^") and ste ep (a = ~0-8) sp ectrum ra- 
dio sources given bylDunlop fc Peacock l)l990l taking the 
MEAN-z data from their appendix C) from 



dN{z) _ dV{z) 



dz 



Plow{z) = X{z) 



X / $,(F,z)dF, 
S \ r 2.7 GHz 



(6) 



(1 + ^) 



l-a 



5.2. The spatial correlation function 

Given the amplitudes of u;(^) determined in Sect.01we can 
use the cosmological Limber equation to estimate the spa- 
tial correlation length, Tq, by deprojecting w{9) into the 
spatial correlation function, ^(r) using the redshift distri- 
bution and cosmology (e.g. Peebles 1980, §56). We con- 
sider two cosmological models: a flat, vacuum dominated, 
low-density Universe (ACDM; 17m = 0.3, 17a = 0.7), and 
an Einstein-de Sitter model Universe (rCDM; ^Im = 1-0, 

= 0). We use Hq = lOO/i km s^^ Mpc'^. 

We assume an epoch dependent power-law spatial cor- 
relation function of the form 



(1 + ^) 



-(3+0 



(7) 



where Vp is the proper distance, tq is the spatial correla- 
tion length'^ at z = 0, and e parameterizes the redshift 
evolution of the clustering. To express ^{rp, z) in terms of 
comoving coordinates rc — rp x {I + z) , we write: 



e(r„z)=(i^) (l + z)^-(3+.)^ 



which can be written as 



, ro{z) = ro(l + z)^ 



(8) 



(9) 



where ro(z) is the (comoving) correlation length measured 
at z. In a flat model Universe, the cosmological Limber 
equation can be expressed as follows (see e.g. IPeeblesI 



Note that the spatial correlation length is not a physical 
lengthscale in the space distribution of galaxies. It is just de- 
fined as that length at which ^(r) is unity (i.e. the chance of 
finding a galaxy at the distance ro from another galaxy is twice 
the Poissonian chance). 
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Fig. 13. Dashed lines show the redshift distributions 
for Si, 4 > 10 mJy, computed from the free-form mod- 
els 1 — 4, the pure luminosity evolution model (PLE) 
and the luminosity /dens ity evolution model [LDE) of 
IPunlon fc PeacocM ljl99f]l) (see text for details) . The av- 
erage of the six different models is indicated by the solid 
curve. 



Il98(il : 




(10) 



with 



Q{z) = [{i + zf + nil-i]''-\ 

, . I r dz 

x[z) - 



(11) 




and using the approximation that angles are small 
1). We calculate N{z) for each subsample. 

The evolution parameter e can represent a variety of 
clusteri ng models. Three iniportant cases are the follow- 
ing fsee lPhillipDS et alJl978tlKund ic'l997^. (1) The stable 
clustering model (e = 0): if galaxy clustering is gravitation- 
ally bound at small scales, then clusters have fixed physical 
sizes (i.e. they will neither contract nor expand) and will 
have a correlation function that decreases with redshift as 
(1-t-z)^"'^'^. (2) The comoving clustering model (e — 7 — 3): 
galaxies and clusters expand with the Universe, so their 
correlation function remains unchanged in comoving co- 
ordinates. This case applies well to a low density Universe 
where there is not enough gravitational pull to counterbal- 
ance expansion, and implies that structures have formed 



0.01 



0.008 



S > 100 mJy 




2 3 
Redshift 



Fig. 14. The redshift distributions for 5*1.4 > 100 mJy. 
See the caption of Fig. ^] for details. 



very early. (3) The linear growth model (e = 7 — 1): clus- 
tering grows as expected under linear perturbation theory. 

Studies of the spatial clustering properties of radio- 
quiet quasars indicate that the clustering history of active 
galaxies, unlike that of normal g alaxies, i s best character- 
ized using a negative value for e. lKundid l)l997|) measured 
the high-redshift quasar-quasar correlation function from 
the Palomar Transit Grism Survey, and found no evidence 
for a decrease in the correlation amplitude of quasars with 
redshift. Moreover, he found that £,qq{z > 2)/£^qq{z < 2) ~ 
1.8, suggesting a.n even highe r amplitude at higher red- 
shifts. Similarly. ICroom et al.l (f2001 'l find almost no evo- 
lution in clustering strength for quasars taken from the 
2dF QSO Redshift Survey out to z ~ 2.5. Therefore, we 
opt for evolution model 2 (i.e. constant clustering in co- 
moving coordinates), which implies e = —1.2 for 7 = 1.8. 
In table 01 we list the results obtained using this model for 
the two different cosmological models. For comparison, we 
also indicate the results using the stable clustering model 
(e = 0). For e = the present-day correlation length is 
1.4 times higher than for e = — 1.2 in both cosmologies. 
However, given the strong peak in the redshift distribution 
at 2: ~ 1, we are effectively measuring clustering at z ^ 1. 
Calculating ro{z 1) in the case of stable clustering using 
Eq- m yields a value that is only ~ 1.1 times lower than 
ro{z ~ 1) = To in the case of e = —1.2. Therefore, the 
value of ro{z ~ 1) is relatively independent of the exact 
value of e. The results for the e = —1.2 (ACDM) case are 
presented in Fig. 1151 We find an approximately constant 
spatial correlation length of ~ 6.0 h^^ Mpc from 3-40 
mJy, compared to ~ 14 Mpc at 200 mJy. 

As we have shown, the possibility that the observed 
flux-dependency of the clustering is just an effect of pro- 
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Table 4. Present-day spatial correlation lengths and la errors derived from the galaxy angular correlation func- 
tion (7 — 1.8) of the NVSS as a function of flux density limit. Listed are the results found using two different 
cosmological models and two different values for the evolution parameter e (see text for details). 





e — 


-1.2 


e 


= 


^low 


rCDM 


ACDM 


rCDM 


ACDM 




ro {h~^ Mpc) 


ro {h~^ Mpc) 


ro {h"^ Mpc) 


ro {h~^ Mpc) 


3 mjy 


3.2 ±0.2 


4.5 ±0.3 


4.8 ±0.3 


6.3 ±0.5 


5 mJy 


3.7 ±0.2 


5.2 ±0.3 


5.6 ±0.3 


7.5 ±0.4 


7 mJy 


3.5 ±0.5 


4.9 ±0.6 


5.3 ±0.7 


7.2 ±0.9 


10 mJy 


3.7 ±0.4 


5.3 ±0.6 


5.7 ±0.7 


7.8 ±0.9 


20 mJy 


3.5 ±0.4 


5.0 ±0.6 


5.5 ±0.7 


7.5 ± 1.0 


30 mJy 


3.7 ±0.4 


5.3 ±0.6 


5.8 ±0.7 


7.9 ±0.9 


40 mJy 


4.3 ± 1.0 


6.1 ± 1.4 


6.7 ± 1.6 


9.1 ± 2.2 


50 mJy 


5.3 ± 1.4 


7.5 ±2.0 


8.2 ± 2.2 


11.2 ±2.9 


60 mJy 


5.0 ±0.7 


7.0 ±1.0 


7.7 ± 1.1 


10.4 ± 1.4 


80 mJy 


5.3 ± 1.2 


7.4 ± 1.7 


8.1 ± 1.8 


10.9 ±2.5 


100 mJy 


5.7 ± 1.1 


8.0 ±1.5 


8.7 ± 1.7 


11.6 ±2.2 


200 mJy 


10.6 ±1.8 


14 ±3.0 


15.4 ±2.6 


19.8 ±3.4 
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Fig. 15. Spatial correlation lengths and la errors de- 
rived from the cosmological w{9) of the NVSS, assum- 
ing an evolution parameter e = —1.2, and the ACDM 
model Universe. The dotted line indicates the flux den- 
sity limit at which FRI- and FRII-type radio sources con- 
tribute roughly equally to 1.4 GHz radio source counts. 
The dashed line indicates the flux density limit above 
which the contribution of FRIIs is > 75%. The top axis 
indicates the effective radio luminosity as a function of 
flux density limit. 



radio sources in the WENSS and GB6 surveys. They con- 
cluded that the clustering signal could change as a func- 
tion of flux density if relatively low and high power radio 
galaxies represent different spatial structures at a similar 
epoch [z ^ 1). T aking the predict e d pop ulation mix of ra- 
dio sources from lWall fc JacksoiJ 1^93), we find that for 
Si. 4 > 10 mJy the fractions of FRIs and FRIIs are about 
equal. However, for 5i.4 > 100 mJy the fraction of FRIIs 
is more than ~ 75%. Given the fractional changes of the 
source populations with flux density limit, the clustering 
amplitudes measured are very well matched by a scenario 
in which the clustering oj powerful radio sources (mostly 
FRII) and average power radio sources (FRI/FRII) are 
intrinsically different, with FRIIs being more strongly clus- 
tered at z ~ 1 than the radio galaxy population on average. 

As pointed out bv iRengelinkI ()l998() and 
iRengelink fc Rottgerind ((1993) the large difference 
in observing frequencies and sensitivities of WENSS and 
GB6 (the limiting 1.4 GHz flux densities probed by these 
surveys correspond to 10 mJy for WENSS and 70 mJy 
for Greenbank, respectively) only allowed them to make a 
comparison between the results, whereas the detection of 
the inferred flux-dependency of rg within a single survey 
would be highly desirable. Our analysis of the clustering 
in the single large-area, intermediate-frequency NVSS 
survey is in agreement with their conclusions. 



jection can be ruled out, since the shape of the redshift 
distribution is relatively constant with flux over several or- 
ders of magnitude (at least above ~ 1 mJy). This automat- 
ically implies that the average radio power of the subsam- 
ples increases with flux density (indicated by the top axis 
of Fig. I15II. An alternative e xpla nation was therefore sug- 
gested bv IRengelinkI l|l99 8h and ' Rengelink fc R ottgerina 
1)1999(1 based on their measurements of the clustering of 



6. Discussion 

6.1. Clustering measurements from literature 

We start this section by making a survey of other cluster- 
ing measurements from literature. However, readers may 
wish to skip directly to Sect. l6.2l for a discussion on these 
measurements and the results presented in this paper in 
their cosmological context. 
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In order to compare results from different studies, all 
values taken from literature were converted assuming a 
fixed slope 7 = 1.8 by setting ro,i.8 = {ro^^^^^-^. All cor- 
relation lengths are expressed in comoving units, and we 
have transformed all valu es to a ACDM cosmology (see 
iMagliocchetti et al.ll200Cl(l . Please note that the list given 
below is not complete, and the reader is kindly invited to 
consult the individual papers and the references therein 
for further information. 



ro = 1.2 ± 0.2 Mpc, a t z = 0.34, z = 0.62, z = 0.97, 
and z — 1.39, respectivelv. ICarlberg et a l. ( 2000) present 
measurements on a sample of L ~ galaxies up to 
z « 0.6 and find a much milder decline from tq = 5.1 ±0.1 
/i-i Mpc at z = 0.10 to ro = 4.2±0.4 h-^ Mpc at z = 0.59. 

Clustering of the local population of I RAS-selected 
galaxies is best fit by ro = 3.4±0.2 h^^ Mpc jFisher et alJ 
I1994I . 



6.1.1. Clusters 

Estimates o f the correlati on length of rich Abell clusters 
are giv en bv lBahcall fc Soncira (1983) and£ostmai^et_al 
' 1992I) who found ro = 24 ± 9 /i"! Mpc. iLahav et al 



1989|) found ro = 21 ± 7 /i ^ Mpc fr om an all-sky s ample 



of t he brightest X-ray clusters, and iDalton et all l)l994 
and ICroft et~al] l|l997l) found ro = 19 ± 5 /i"^ Mpc and 
ro = 16 ± 4 Mpc, respectively, for clusters selected 
from t he APM Galaxy Survey. Recently, .Gonzalez et all 
ll2002l) measured the correlation length of distant clusters 
in the Las Campanas Distant Cluster Survey and found a 
correlation length of 24.8 ± 4.5 Mpc at z = 0.42. 

Different studies may have sampled clusters of differ- 
ent degrees of richness, which can account for most of the 
scatter in the reported values. In general, however, all re- 
sults are consistent with clusters being the most strongly 
clustered objects known in the Universe. 

6.1.2. Optically-selected ordinary galaxies and IRAS 
galaxies 

Bright early-type galaxies are found to h ave a strongly 
cluste red distribution in the local Universe. I Willmer et al.l 
lll99j^ find r n = 6.8 ± 0.4 fe-^ M pc for local L > el- 
lipticals, and lGuzzo et a D ll1997^ measure a considerably 
higher ro = 11.4±1.3/i~^ Mpc for a sample of similar 
galaxies. Although these results are only consistent with 
each other at the 3a level, the latter sample contains a 
higher fraction of local clusters, presumably responsible 
for boosting the ro. The dependence of galaxy cluster- 
ing on luminosity and spectral type has been studied us- 
ing the ong oing 2 degree Field G alaxy Redshift Survey 
(2dFGRS). iNorberg et all (|2002') find ro = 11.8 ± 1.6 
Mpc for the brightest early- type galaxies in the 2dF- 
GRS. Moreover, they find a strong dependence of cluster- 
ing strength on luminosity, with the amplitude increasing 
by a factor of 2.5 between and 4L,. The ordinary 
population of galaxies has been found to be less strongly 
clustered t han the population consisting of local (bright) 
ellipticals: iLovedav et aP l(l99 5) find ro = 4.7 ± 0.2 h-^ 
Mpc from the APM survey. At higher redshifts, the clus- 
tering strength in a sample of faint /-C-selected galaxies 
with minimum rest-frame luminosities of = —23.5, or 
about O. SLa., is found to be fa irly rapidly declining with 
redshift: ICarlberg et aTl ll1997^ find ro = 3.3 ± 0.1 
Mpc, ro = 2.3 ±0.2 h'^ Mpc, ro = 1.6 ±0.2 Mpc, and 



6.1.3. Extremely red objects (EROs) 

Several recent studies indicate that the comoving cor- 
relation length of early-type galaxies undergoes little or 
no evolution from < z < 1. Evidence for this is pro- 
vided by the clustering of extremely red objects, a pop- 
ulation of galaxies having very red optical to infrared 
colors {R — Kg > 5). These red colors are consistent 
with them being either old, passively evolving elliptical 
galaxies, or strongly dust-enshrouded starburst galaxies at 
z ~ 1 — 1.5. Indeed, further observations have confirmed 
that both classes are present in the ERO populatiori (e.g 



that both classes are present m the bKU populatiori (e.g. 
Dunlop et al.lh 996: Ci matti et al.lll99Sl IPev et a,lJ h99i 
Liu et alJ I2OOQI . .Daddi et alJ ()200ll) have recently em- 



barked on a study of the spatial clustering of a large sam- 
ple of i > EROs at z ~ 1, and found a large corre- 
lation length of ro = 12 ± 3 Mpc. In ICimatti et alj 
1^2002) the results are presented involving the EROs that 
were identified in a large flux limited redshift survey of 
~ 500 galaxies with K < 20. The derived fraction of 
early-type EROs from that sample is 50±20%, while there 
is an increasing contribution of dusty star-fo r ming EROs 
at faint mag nitudes. Therefore, iDaddi et all l|2002 ') have 
attempted to analyse separately the spatial clustering of 
EROs from both categories by studying the frequency of 
close pairs. They find that the comoving correlation length 
of the dust-enshrouded starbursts is constrained to be less 
than ro = 2.5 Mpc, while the old EROs are clustered 
with 5.5 ^ rg < 16 h"^ Mpc. This is consis tent with 
the value reported earlier in iDaddi et al.l ((200 ll) , which is 
still valid as a lower limit for the clustering of early-type 
EROs based on the argument that the much less clus- 
tered dusty star-forming EROs only dilute the clustering 
signal coming from the ellipticals in this sample (see also 
iRoche et al.ll2n02j) . Furthermore. iMcCarthv et all ((200lh 
have identified a large sample of such faint red galaxies as 
being consistent with mildly evolved early-type galaxies 
at z ~ 1.2. They find a clustering strength of ro = 9.5 ± 1 
h-^ Mpc. 



6.1.4. Radio galaxies 

The results on the spatial clustering of radio sources at 
z ~ 1 presented in this paper indicate that ro depends on 
radio luminosity in such a way that very luminous (FRII) 
radio galaxies cluster more strongly than the total popu- 
lation of radio galaxies (both FRI and FRII) on average, 
reminiscent of a similar luminosity trend found for samples 
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of optically-selected galaxies. We roughly construct two 
radio luminosity bins from our measurements by compar- 
ing the tq found for 3-40 mJy to the Tq found for the 200 
mJy subsample. We find tq ~ 6±1 Mpc~^ for the rel- 
atively low power bin (Pi. 4 ~ io24-25 -y^ Hz^^ sr^^), and 
ro ~ 14± 3 h-^ Mpc~i for the high power bin (P > 10^*^ 
W Hz-i sr-i). 

6.1.5. Optically-selected quasars 

ICroom et all l|200l|) have determined the correlation 
length of quasars (QSOs) using 10,558 quasars taken from 
the 2dF QSO Redshift Survey. They find that QSO clus- 
tering appears to vary little with redshift, with ro — 
4.9 ± 0.8 /i-i Mpc at z = 0.69, tq = 2.9 ± 0.8 h-^ Mpc at 
z = 1.16, ro = 4.2±0.7 /i"! Mpc at z = 1.53, tq = 5.3±0.9 
h-^ Mpc at z = 1.89, and ro = 5.8 ± 1.2 h''^ Mpc at 
2 = 2.36. 



6.1.6. Lyman-break galaxies 

Lyman-brcak galaxies (LBGs) are found to be associated 
with star- forming galaxies at z ~ 3, with co moving cor- 
relati on lengths of ro = 3.3 ± 0.3 Mpc llAdelberger 
2000t ). and ro = 3.6 ±1.2 h''^ Mpc l)Porciani fc Giavaliscc 
2002). Ouchi ct al. (2001) find ro = 2.7±0.6 h'^ Mpc for 
a sample of LBGs at z ~ 4. 



6.2. Clustering evolution 

6.2.1. The clustering of massive ellipticals at z ~ 1 

In Fig. 1161 we present an overview of the evolution of 
galaxy clustering, as it follows from the broad variety of 
observational results summarized above. The ro that we 
measure for the brightest radio sources at z 1 is com- 
parable to the ro measured for bright ellipticals locally, 
and ^ 2 X higher than the ro measured for relatively faint 
radio sources and quasars, suggesting that they are con- 
siderably more biased and probably probe spatial struc- 
tures associated with strongly clustered, massive objects. 
This does not come totally unexpectedly, as there is a 
range of observat ional evidence in support of this result. 
iBest et al] l|l998|) found that powerful 3CR radio galax- 
ies are mostly associated with massive galaxies at z ~ 1, 
and at high (z ~ 1) and very high (z > 2) redshifts the 
most luminous (i.e. FRII-type) radio sources are found 
in very dense environments associated with forming clus- 
ters. This is based on for example the presence of large 
X-ray halos llCra.wford fc Fabian Il996l). excesses of com- 
panion g alaxies ('Mc Carthvlll988nR,ottgering eraP Il996l: 
iNakata et al. 20011 . an d excesses of Lya emitters around 
powe rful radio sources l)Kurk et a,].ll2nnfil:IVenema,ns et a,l.l 
12002(1 ■ Furthermore, most very high redshift radio galaxies 
(z > 2) are surrounded by g i ant halos of emission lin e gas 
(e.g. lR5ttgering et alJ 119991 iDe Breuck et al.ll2000|) . and 
some have very clumpy morphologies suggestive of massive 
forming systems ('e.g. IPentericci et al.iil999L .2000,1 . Using 



HST/NICMOS observations. IPentericci et alJ (|2001^ have 
found a number of radio galaxies at z ~ 2 having mor- 
phologies that are represented well by a de Vaucouleurs 
profile, consistent with them being elliptical galaxies or 
proto-galaxy bul ges. 

As argued bv lBest et al.1 ()l999(l . powerful radio sources 
must rely on (i) a plentiful supply of gas to fuel a super- 
massive blackhole that can drive the AGN activity, and 
(ii) a dense surrounding medium able to contain the radio 
lobes. These environments are indeed expected to be found 
in the gas-rich galaxy clusters at high redshift, addition- 
ally supporting the conclusion that high redshift FRIIs are 
associated with strongly clustered, massive objects. One 
may argue that this conclusion somewhat contradicts the 
fact that low redshift FRIIs are primarily found to be sit- 
uated in small, isolat ed galaxy groups, and not in the cen- 
ters o f large clusters (|Butcher fc Oemleill978tlHill fc Lillvl 
I1991I) . This, however, can easily be explained by consider- 
ing that the local analogs of the gas-rich cluster environ- 
ments that are suitable for producing powerful FRIIs at 
high redshifts, are found in relatively small galaxy groups, 
and not in the g as-depleted centers of local rich clusters 
( Rengefinklll998|i . 

Interestingly, we find that both EROs and powerful ra- 
dio gala xies are strong l y clus tered with ro > 10 Mpc 
at z ~ l. lWillott et"all|200lh suggested that high-redshift 
radio galaxies and EROs could be identical galaxies seen at 
different stages of their evolution, based on their findings 
of ERO-like host galaxies for a number of radio galaxies 
from the 7C Redshift Survey. This, of course, would be 
highly consistent with the belief that both radio galaxies 
and EROs may be the progenitors of local bright ellip- 
ticals. They conclude that the density of radio sources 
with minimum radio luminosities of logio P151 = 24 W 
Hz^^ sr^^ is consistent with a model in which all EROs go 
through a relatively short period of AGN activity, forming 
a radio galaxy somewhere between z = 2 and z = 1. 

However, if all EROs are radio galaxies at some stage, 
their highly clustered spatial distribution should be re- 
flected in the spatial distribution of the radio galaxies. Fig. 
Elshows that the clustering of EROs and radio galaxies is 
consistent only for those galaxies with radio luminosities 
of logio Pi4oo 26 W Hz^-^ sr^^. The surface density of 
such radio sources in the redshift range 1 < z < 2 in the 
NVSS is ~ 2 X 10~* arcmin"^, while the surface density of 

llnaddi et ahlEf 

ticals among EROs to be 



2001). If we take the fraction of old elliT 



i d enip - 

70% (ICima,tti et a1.ll200^ . 

then only ^ 0.06% of these EROs are currently observed 
in their radio-loud phase. However, because the typically 
assumed AGN lifetimes are short compared to the cosmo- 
logical time-scale from z = 2 and z = 1 (^^=2-1 — 3.5 
Gyr for f^Af = 0.3, f^A = 0.7), the number of EROs that 
could undergo a radio- loud phase is ~ 2 — 20% (assuming 
tAGN — 10''"® yr.). These fractions can be increased sig- 
nificantly if, for example, we select EROs that are much 
redder: the density of EROs having R — Kg > 6 is a fac- 
tor of ~ 10 lower compared to i? — if^ > 5 ( Daddi et alJ 
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Fig. 16. The redshift evolution of galaxy clustering in a ACDM Universe. See the text for references to data taken 
from literature. Lines represent the following models: (i) stable clustering (e = 0) normalized to tq of local ellip- 
ticals and clusters {dotted lines), (ii) Hnear clusteri ng (e = 1) normaliz ed to 5 h~^ Mpc (dot-dashed lin e ), (iii) 
clustering of the daA matter (tMck soM Une, from I .T^nkins et alJ l|l998t) , see also iMoustakas &: Somervilld l|2002li 
for a use ful p arameterization), (iv) galaxy conservation model normalized to rg of local ellipticals {thin solid line, 
see fFrvl |l99a) ). (v) hierarchical r nodel for clustering evolu tion of early- type galaxies normalized to rg of local el- 
lipticals {thick dashed line, from iKauffmann et all l|l999lJl ). and (vi) clusterin g evolution as a function of dark 
matter halo masses with ^Amin = 10^^"^'' M© {thin dashed lines, from [Matarr ese et alJ l|l997l) '). A nice represen- 
tation of this figure showing actual images of the various objects rather than symbols can be found at our website: 
http://www.strw.leidenuniv.nl/~overzier/rO.html, 



l200ll) . giving 14 — 140%. It may be clear from the above 
that the unification of EROs and radio galaxies, although 
tempting, relies on a number of issues that have not yet 
been resolved. Further study of the luminosities, colors 
and morphologies of radio galaxy hosts, as well as the 
cluster environments of EROs may be expected to pro- 
vide important clues for constraining this scenario. 

6.2.2. Comparison with theoretical predictions 

Linear (e ~ 1, dot-dashed line) or stable (e = 0, dot- 
ted line) clustering evolution models have been found to 
best fit the measurements of ordinary, optically-selected 



galaxies at z < 1 (e.g. ICarlberg et~el\ Il997l l2000t 
iMcCracken et al. l200l[ and references therein). However, 
as Fig. El shows, these models do not provide a good de- 
scription for the evolution of massive early type galax- 
ies as inferred from the meassurements of local bright 
ellipticals and FRII radio galaxies and EROs at z ~ 1. 
Adjusting these models to the measurements would ei- 
ther require z ~ 1 massive ellipticals to have a correlation 
length around 6 — 7 ft." ^ Mpc, or local bright ellipticals 
to have a correlation length of the order of that of local 
clusters, far greater than observed. For these galaxies, the 
current measurements require a model that predicts rela- 
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tively constant clustering in comoving coordinates, i.e. a 
negative value of e « — 1 in the simple e-model. 

Although the parameterization of clustering evolution 
by means of the e-model is useful for characterizing the 
measurements as a function of redshift, it does not pro- 
vide good physical insight into evo lution governed by the 
clustering of dark matt er halos Csee lGiavalisco et al 
iMcCracken et al.ll200l|) . Galaxy clustering evolution can 
be described more precisely by 



^gai{z,r)=D^{z)b\z)^rn{0,r), 



(12) 



where D(z) is the linear cosmological growth rate (see 
ICarroll et alJll99^ . b{z) the evolution of the bias, and 
Cm(0, r) the correlation function of the underlying matter 
distribution at z = 0. Since b{z) is related to the nature of 
the mechanism through which the galaxies were formed, 
measurements of ^gai{z, r) can be used to constrain struc- 
ture formation models. 

In the galaxy conservation model, objects are formed 
by means of monolithic collapse at arbitrarily high red- 
shift, and their clustering evolution is described solely 
by th e cosmological growth of density perturbations (^Frvl 
Il99fil) . In this model, bias evolves as 



b{z) = l + {b^-l)/D{z), 



(13) 



where 6o = (crs.gai/crs.m) and cg is the rms fluctuation 
amplitude inside a sphere of 8 Mpc radius. Taking 
^"0,771(0) = 5 Mpc for the present-day correlation 
length of the d ark matter from the GIF /VIRGO iV-body 
simulations of I.Tenkins et alJ l(l998l) (thick solid line in 
Fig. Unj and fo,gai(0) = 8.7 h~'^ Mpc for eUipticals, 
we find cg^m — 0.9 and 

(^s>,gai — 1.5 corresponding to 
60 ~ 1.65. This model is indicated in Fig. ^| (thin solid 
line). Analogous to the above arguments against simple 
stable or linear clustering, extrapolating the clustering of 
local ellipticals to z ~ 1 in the galaxy conservation model 
does not fit the observed extreme clustering of EROs and 
powerful radio galaxies. On the other hand, this scenario 
shows good agreement with the r^ ^ (S h"^ Mpc measured 
for lower luminosity radio sources and QSOs at z ~ 1. 

Crucial to t he picture tha t is developing may be the 
recent results nf lWilsonl ll2nn,Sh. who studied the clustering 
of {V — /)-selected L, early-type galaxies in the redshift 
range 0.2 < z < 0.9. This author found that these galaxies 
cluster slightly more strongly compared to the field, with a 
best-fitting e-model of e = and ro = 5.25±0.28 h'^ Mpc. 
This is in agreement with the correlat i on len gth of local 
early- types in the 2dFGRS. 'Wilson' (2003) remarks that 
this measurement is inconstent with the large found 
for EROs, which are also believed to be L ~ _L» early- 
type galaxies. The value of rg for EROs and radio galax- 
ies could be spuriously high due to uncertainties in their 
redshift distributions which is not included in the quoted 
errors, although the selection functions of both EROs and 
powerful radio gal axies are considered to be understood 
relatively well (e.g. iDunlop fc Peacocklll990t iDaddi et alJ 



because they correspond to a population of massive, bright 
cluster galaxies in the process of formation. If FRII radio 
galaxies and EROs are indeed the distant analogs of lo- 
cal L ~ L* early-types, they are becoming considerably 
more biased tracers of the underlying galaxy distribution 
with redshift, while this galaxy distribution itself prob- 
ably traces the dark matter distribution with relatively 
constant bias. Interestingly, (semi-) analytic models and 
iV-body simulations are able to explain this bias evolu- 
tion and the large inferred rp at z '--^ 1 of massive ellipti- 
cals, if the assumption that galaxies are conserved quanti- 
ties (i.e. closed-box s ystems) is relajced. These hierarchical 
merg ing models (e.g. IMo fc Whitelll996t IMatarrese et all 



1997t iMoscardini et all Il998l: iKauffmann et al l Il999hl: 
Moustaka,s fc Somervil]el2002HMo fc Whitel20oTa,nd ref- 



erences therein) prescribe that for certain types of objects 
bias can grow stronger with redshift than the growth of 
perturbations, resulting in a rp that is constant or even 
increasing with redshift. 

In the (transient) model of lMatarrese et al] ( 199?!) it is 
assumed that the mass of the dark matter halo also deter- 
mines the physical paramet ers of the galaxy tha t it con- 
tains. Based on the work of iMo fc Whit j lll99fil) and the 
forma lism of IPress fc Schechteil l)l974|) . IMatarrese et alJ 
derive that the bias in such a model evolves as 



b{z) = 1 - 1/5, + [60 - (1 - \l^c)\lD{z)^ 



(14) 



where be = 1.686 is the critic al linear o verdensity for 
spherical collapse (but see also lLiii3 ll992l) . The param- 
eters &o and /3 depend on the minimum mass of the halo, 
and we have used the COBE-no rmalized (ACDM) values 
for A^min = 10^2-14^^ gj.^gj^ bv lMoscardini et all l|l998l) 
to plot this model in Fig. ^|(thin dashed lines). We find 
that the Almin = VS^'^Mq model is able to fit the mea- 



surements at both 



1 and z = 0. Likewise, the model 



l200lUMcCarthv 
dio galaxies at z 



il.l l200lD . Alternatively, EROs and ra- 
1 may be much more strongly clustered 



with A^min = 10^^ Af© has been f ound to fit the sp atial 
clustering of QSOs relatively well (iCrooni et alJl200ll). al- 
though several seri ous caveats exist fsee iRengelinkHlQQSt 
ICroom et all 120011) . Most importantly, the assumption 
that there always exists a simple relationship between the 
mass of the dark matter halo and the property by which 
a galaxy is selected may not be valid. 

In Fig. ^1 we have also indicated the predicted evo- 
lution of the clustering of early-typ e galaxies (thick 
dashed line) from t he ACDM-models of IKauffmann et alJ 
(tl999a,b) (sec also ISomerville et al.ll200lh . normalised to 
ro found for local el lipticals. An impo r tant fe ature of the 
models presented in *Kauff mann et aD 1 1999bl) is that one 
naturally expects a dip in rg between z = and z « 1, 
if structure is probed by galaxies of intermediate lumi- 
nosities residing in haloes of masses 10"'^^~^^Mq that have 
formed early and are unbiased tracers of the overall mass 
distribution. However, these simulations also show that 
this dip is very sensitive to sample selection criteria: mas- 
sive early-type galaxies exhibit no dip in clustering be- 
tween z — and z « 1, because they occur in rare, 
very masssive haloes of 10^^"^'* M© which are strongly bi- 
ased locally, and which become even stronger biased with 
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redshift. The agreement of this model wit h the results 
presen ted i n this paper and the results of iDaddi et alJ 
l|200l[) and IMcCarthv eraP (|20o3) is striking. Although 
promising, some discrepancies betwe en the mod e l and 
the o bservations remain. For instance. iDaddi et alJ l|200lL 
find strong disagreement between the model and 
the high observed space density of EROs, seemingly con- 
sistent with the purely passive evolution of local ellipti- 
cals. Furthermore, current merging models generally pre- 
dict that these galaxies should have experienced recent 
star-formation activity, while this is not observed. It may 
become possible to still reconcile the observations with 
the ACDM merging models if, for exar nple, the mer ging is 
accompanied by little star- formation ijPaddi et alj EoOll. 
Also, the EROs are found to have relatively old stellar 
populations of > 3 Gyr that sho w no indications of recent 
forma tion processes. However, iMoustakas fc Somervillj 
l|200 2^ point out that the relatively old ages of their stel- 
lar populations do not automatically imply similar ages 
for the host galaxies. 

Despite the success of current hierarchical models in 
predicting the evolution of bias for these massive galax- 
ies, we would like to point out that galaxy conservation 
or linear/stable clustering evolution could still be able to 
explain the measurements if EROs and/or powerful radio 
galaxies are solely found in rich Abell-type clusters with 
(present-day) tq ~ 15 — 25 Mpc. As we have shown 
there is substantial evidence that this may be the case 
for, at least, the powerful radio galaxies, and future data 
may show whether this also holds for (a subset of) the 
population of EROs. 

At the highest redshifts, clustering of LBGs at 3 < 
z < 4 indicate that these objects can be connected to lo- 
cal ellipticals in a galaxy conservation scenario. However, 
it is now believed that LBGs probably occupy much less 
massive halos of 1O^^~^^M0 than those that contain lo- 
cal massive galaxies, suggesting that if these objects are 
to be the progenitors of local ellipticals, the y must have 
accumulated a considerable amount of mass I Adelbergerl 
E00(ilMousta,ka,s SomerviliaEoO^. 



6.2.3. Clustering and the occurrence of AGN at high z 

Fig. 1161 suggests that the clustering evolution of active 
galaxies in general is considerably different from that of 
ordinary galaxies. Albeit at a lower amplitude, the cluster- 
ing of QSOs also shows a trend of constant or slightly in- 
creasing amplitude with redshift, very similar to the trend 
that we derive for the clustering of the most massive el- 
lipticals. According to the standard paradigm, AGN are 
powered by the accretion o f matter onto a (super-)massive 
blackhole fe.g. lR,ee'3ll984l) . This fuelling mechanism may 
very well be associated with the injection and accretion 
of gas during major merging events, and thus, the occur- 
rence of AGN seems to be logically linked to the hierarchi- 
cal scenari os for structure formation. Recently, in a series 
of papers IjHaehnelt et al.i.1998; .Haehnelt fc Kaiiffmanni 



12000^'Ka uffmann fc Haehneltll200(l l2002|) the simulations 
of Kauff mann et alJ l|l999bf) were extended to a unified 
model for the evolution of both galaxies and quasars. In 
their model, elliptical galaxies, supermassive black holes 
and starbursts are formed during major merging events, 
in which a fraction of the available gas is used to trig- 
ger quasar activity by accretion for about 10^ years, and 
the remaining gas is converted into stars in a single short 
burst. This model succesfuUy reproduces the evolution of 
cold gas that is derived from observations of damped Lya 
systems, the luminosity functions and clustering proper- 
ties of QSOs from the 2dF QSO survey, and the relation 
between bulge velocity dispersion and black hole mass 
that has been found in den iographic studies of black holes 
in nearby galaxies ( e.g. iKormend v fc Ric hstonel Il995t 

iMagorrian et al.lll998l: iGebhardt et al 200m . 

Although it has yet remained unknown exactly what 
processes cause the physical differences between radio- 
quiet and radio-loud AGN, recent results indicate that 
the hosts of all powerful AGN (both radio-loud and radio- 
quiet) are almost excl usively L > L^, ellipticals (see 
iDunlop fc McLurel2003l and references therein). However, 
the same studies also indicate that while radio-quiet AGN 
hosts can have black holes with masses of 1O^~^°M0, 
the radio-loud sources are cleanly confined to black hole 
masses M^h Jc^ 5 x 1O*M0. Furthermore, in the regime 
of extreme radio luminosities that lie well beyond the 
FRI/FRII luminosity-break, the power needed can only 
be achieved by blackholes with Afb/i > IO^Mq, requiring 
host masses of > IO^^A/q that imply L > L* luminosities 
l)DunloD fc McLure-2003ft . This may explain why the most 
powerful NVSS sources are extremely clustered compared 
to the, on average, less massive hosts of QSOs. This is 
supported by the fact that the radio sources in our lower 
radio luminosity bin have a correlation length similar to 
that of QSOs at z ^ 1, while both populations are still 
clustered more strongly compared to the field at z ~ 1. 
We conclude that the masses of the haloes, host galax- 
ies, and black holes that are probed by the most powerful 
radio sources are among the most massive objects in the 
Universe, possibly formed through massive mergers in hi- 
erarchical fashion. 

7. Summary 

The main conclusions that can be drawn from our 
analysis are the following: 

• Below ^ 6' w{d) is dominated by the size distribution 
of multi-component radio sources. A simple model of the 
physical size distribution of FRII radio galaxies is able 
to explain the observed enhancement of the cosmological 
clustering signal. 

• The amplitude of the angular two-point correlation 
function of radio sources increases with increasing radio 
flux, corresponding to a similar increase in with 
increasing average radio power of the samples. This sug- 
gests that powerful FRII radio galaxies are intrinsically 
more strongly clustered than the average population 
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of radio galaxies at 2; ~ 1. This is consistent with the 
extremely rich environments in which high redshift FRIIs 

are generally found. 

• The correlation lengths of powerful radio galaxies 
and EROs are of comparable magnitude and both are 
associated with massive ellipticals at z ^ 1. This suggests 
that we could be looking at identical objects at different 
stages of their evolution, implying that AGN activity is 
an important phase in the evolution of massive galaxies 
in general. 

• The evolution that we infer for the clustering of massive 

ellipticals between z ^ 1 and z ^ is in agreement 
with predictions from hierarchical models for structure 
formation, because they can account for the observed lack 
of evolution in rg. However, the large correlation length 
of powerful radio galaxies at 2; ~ 1 is also consistent 
with galaxy conservation models if they are primarily 
associated with rich, Abell-type clusters. 
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